On blowup for Yang-Mills fields 
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We study the development of singularities for the spherically symmetric Yang-Mills equations 
in d + 1 dimensional Minkowski spacetime for d = 4 (the critical dimension) and d = 5 (the 
lowest supercritical dimension). Using combined numerical and analytical methods we show in both 
cases that generic solutions starting with sufficiently large initial data blow up in finite time. The 
mechanism of singularity formation depends on the dimension: in d = 5 the blowup is exactly 
self-similar while in d = 4 the blowup is only approximately self-similar and can be viewed as 
the adiabatic shrinking of the marginally stable static solution. The threshold for blowup and the 
connection with critical phenomena in gravitational collapse (which motivated this research) are 
also briefly discussed. 
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Introduction: The Yang-Mills (YM) equations are the 
basic equations of gauge theories describing the funda- 
mental forces of nature so understanding their solutions 
is an issue of great importance. This is not an easy 
task since, in contrast to Maxwell's equations or the 
Schrodinger equation, the YM equations are nonlinear 
which opens up the possibility that solutions which arc 
initially smooth become singular in future. Actually such 
a spontaneous breakdown of solutions of YM equations 
cannot occur in the physical 3-1-1 dimensional Minkowski 
spacetime as was shown in a classic paper by Eardley 
and Moncrief who proved that solutions starting from 
smooth initial data remain smooth for all future times. A 
natural question is: how does the property of global reg- 
ularity depend on the dimension of the underlying space- 
time, in particular can singularities develop in d + 1 di- 
mensions for d > 3? We hope that our letter is a step 
towards answering this question. As we argue below, 
the problem of singularity formation for YM equations in 
higher dimensions is not only interesting in its own right 
but in addition it sheds some light on our understanding 
of Einstein's equations in the physical dimension. 

Despite intensive research the problem of global regu- 
larity for YM equations in 4 -|- 1 dimensions is entirely 
open [||. A lot of progress has been made to prove lo- 
cal existence for "rough" initial data, yet the attempts 
of proving global regularity by establishing local well- 
posedness in the energy norm fail to achieve the goal by 
"epsilon" (nota bene such a local proof of global exis- 
tence has been obtained in 3 -I- 1 dimensions , thereby 
improving the theorem of Eardley and Moncrief). In this 
letter we report on numerical simulations which in combi- 
nation with analytic results strongly suggest that generic 
solutions with sufficiently large energy do, in fact, blow 
up in finite time. Hence, we believe that the above men- 
tioned epsilon in the optimal local well-posedness result 
is not a technical shortcoming but is indispensable. We 



show that the singularity formation is due to concentra- 
tion of energy and has the form of adiabatic shrinking of 
the marginally stable static solution. 

Higher dimensions, d > 4, appear to be somewhat 
under-explored; the only result we are aware of is the 
proof of existence of self-similar solutions in d = 5, 7, 9 
These solutions provide examples of singularities develop- 
ing from smooth initial data, however nothing was known 
about their genericity and stability. Here we restrict our- 
selves to the d = 5 case because of its connection with 
Einstein's equations. We first show that the example of 
self-similar blowup given in Q is, in fact, generic. Then 
we look at the threshold for singularity formation and 
observe a behaviour similar to the critical behaviour in 
gravitational collapse |^ (with blowup being the ana- 
logue of a black hole), in particular we find a self-similar 
solution with one instability as the critical solution. 

We remark in passing that there are close parallels be- 
tween YM equations in d -I- 1 dimensions and wave maps 
in (d— 2) + l dimensions Indeed, many of the phenom- 
ena described below have been previously observed by us 
for the equivariant wave maps into spheres in two |^] and 
three ^ spatial dimensions. 

Setup: We consider Yang-Mills fields in d -I- 1 dimen- 
sional Minkowski spacetime (in the following Latin and 
Greek indices take the values 1, 2, . . . , d and 0, 1, 2, . . . , d 
respectively). The gauge potential Aa is a one- form with 
values in the Lie algebra g of a compact Lie group G. 
Here we take G — SO{d) so the elements of 5 = so(d) 
can be considered as skew-symmetric d x d matrices and 
the Lie bracket is the usual commutator. In terms of the 
curvature — daAp — dpAa + e[Aa,Ap\ the Yang-Mills 
equations are 
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where e is the gauge coupling constant. It is customary 
to set e = 1 and we shall also do so in the following. 
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However, it is worth remembering that [e^] = M~^L'^~'^ 
(in c = 1 units) ; in particular has the same dimension 
in d = 5 as Newton's constant G in d = 3. 

The YM equations (|l]) are scale invariant: if Aa{x) is 
a solution, so is Aa{x) = Aa{x/ \). The conserved 
energy 
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scales as £{A) = A''^'*£(j4), thus in the PDE terminology 
the YM equations are subcritical for d < 3, critical for 
d = 4, and supercritical for d > 5. It is believed that 
subcritical equations are globally regular because energy 
conservation rules out concentration of solutions on ar- 
bitrarily small scales. In contrast, for supercritical equa- 
tions concentration might be energetically favourable and 
consequently singularities are expected to occur. The 
critical equations provide an interesting borderline case. 
We assume the spherically symmetric ansatz Q 
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where r = \/x^- Then, the YM equations (|l|) reduce 
to the scalar semilinear wave equation for the magnetic 
gauge potential w{t,r) 
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center, reaches a minimal radius, bounces back and then 
disperses to infinity leaving behind an empty space. For 
large amplitudes the ingoing wave keeps concentrating 
near the center and eventually blows up in finite time. As 
the blowup time T is approached we observe the develop- 
ment of a rapidly evolving inner region which is clearly 
separated from an almost frozen outer region. The in- 
ner solution attains a kink-like shape which shrinks in a 
self-similar manner 
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where the profile W depends on the dimension but oth- 
erwise seems universal. The scale X{t) goes to zero as 
t ^ T which signals blowup since the second derivative 
d^w{t,0) = W"{0)X~'^(t) becomes unbounded (there is 
no blowup of the first derivative because W'{0) = 0, as 
we shall see below). 

The subsequent discussion of the details of blowup has 
to be given separately in each dimension. We begin with 
the easier supercritical case. 

d = 5: In this case the scale changes linearly, that is 
X(t) = T — t, as one would expect from dimensional anal- 
ysis. The blowup profile is given by the exact self-similar 
solution of equation (0) 
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where A((j_2) = df- + ^^dr is the radial Laplacian in 
d — 2 dimensions. We solve numerically the initial value 
problem for the above equation in d = 4 and 5. Our sim- 
ulations were performed using finite-difference methods 
combined with adaptive mesh refinement. The latter was 
essential is resolving the structure of singularities devel- 
oping on very small scales. To ensure regularity at the 
center we require that 0) = 1 + O(r^). At the outer 
boundary of the computational grid we impose the out- 
going wave condition. Below we present results for the 
time-symmetric gaussian initial data of the form 

w{0,r)^l-Ar^exp[-(j{r-R)^], wt{0,r)^0, (5) 

with adjustable amplitude A and fixed parameters a = 10 
and R — 2. We have obtained the same qualitative re- 
sults for several other families of initial data so we believe 
that the phenomena described here are generic. 

Results: We begin our description in a unified 
dimension-independent manner; all statements which do 
not cxplicitely involve the dimension apply both to d = 4 
and d ~ b. Since our data are time-symmetric, the initial 
profile splits into ingoing and outgoing waves. The evo- 
lution of the outgoing wave has nothing to do with singu- 
larity formation so we shall ignore it in what follows. The 
behaviour of the ingoing wave depends on the amplitude 
A. For small amplitudes the ingoing wave approaches the 
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FIG. 1. Blowup in d = 5. The plot shows the late time 
evolution of large initial data {A = 0.2). As the blowup pro- 
gresses, the inner solution gradually attains the form of the 
stable self-similar solution W(i{r /{T — t)). The outer solution 
appears frozen on this timescale. 

The solution Wq was proved to exist in and recently 
found in closed form by one of us [|o|. This solution is 
linearly stable ||l0[| (apart from the instability correspond- 
ing to shifting the blowup time) which supports the fact 
that we see it as a generic attractor without tuning any 
parameters of initial data (see Fig. 1). 
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We think that the basic mechanism which is responsi- 
ble for the observed asymptotic self-similarity of blowup 
can be viewed as the convergence to the lowest "en- 
ergy" configuration. To see this, note that rewriting (^) 
in terms of the similarity variable rj and the slow time 
T = — ln(T — t) one can convert the problem of blowup 
into the problem of asymptotic behaviour of solutions 
for T — > oo. The point is that in these variables the wave 
equation contains a damping term which simply reflects 
the presence of an outward flux of energy through the 
past light cone of the singularity. Hence it is natural to 
expect that solutions will tend asymptotically to the least 
"energy" equilibrium state, which is nothing else but Wq. 

It was shown in [^o) that Wq is actually the ground 
state of a countable family of self-similar solutions Wn 
(n = 0, 1, . . .) of equation (^). All n > solutions are 
unstable and therefore not observed in the evolution of 
generic initial data. However, they may show up in the 
evolution of specially prepared initial data. The solu- 
tion Wi with one unstable mode is particularly inter- 
esting since it appears as a transient metastable state 
in the evolution of initial data tuned to the threshold 
for blowup. This indicates that Wi is a critical solution 
whose codimension-one stable manifold separates blowup 
from dispersion (see Fig. 2). 
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FIG. 2. The critical behaviour in d = 5. The rescaled so- 
lution w{t, {T — t)r) is plotted against ln(r) for a sequence of 
intermediate times. Shown (solid and dashed lines) is the 
pair of solutions starting with marginally critical amplitudes 
A = A*±e, where A* = 0.144296087005405. Since e = 10"^^ 
the two solutions are indistinguishable on the first seven 
frames. The convergence to the self-similar solution Wi (dot- 
ted line) is clearly seen in the intermediate asymptotics. The 
last two frames show the solutions departing from the interme- 
diate attractor towards blowup and dispersion, respectively. 

Let A* be the critical amplitude corresponding to the 
threshold. For initial data with amplitudes near A*, in 
the intermediate asymptotics when the solution hangs 
around Wi, the amplitude of the unstable mode about 



Wi is proportional to {A* — A){T — t)~'' where 7 = 5 
(!) is the eigenvalue of the unstable mode This im- 
plies that the time of departure from the intermediate 
attractor, call it t*, scales as T - \A* - Aj^/^. Vari- 
ous scaling laws can be derived from this. For example, 
consider solutions with marginally subcritical amplitudes 
A = A* — e. For such solutions the energy density 
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initially grows at the center, attains a maximum at a 
certain time tmax and then drops to zero. An elementary 
dimensional analysis based on the above scaling predicts 
the power law p{tmax,0)] ^ e"**/^. We have verified this 
prediction in our simulations (with 4% error). 

We point out that the threshold behaviour described 
above shares many features with critical phenomena at 
the threshold for black hole formation in gravitational 
collapse [||. This fact, together with similar results for 
wave maps in 3 -f- 1 dimensions p| , (and other sys- 
tems ) shows that the basic properties of critical col- 
lapse, such as universality, scaling, and self-similarity, 
originally observed for Einstein's equations, actually have 
nothing to do with gravity and seem to be robust prop- 
erties of supercritical evolutionary PDEs. 

d = 4: In this case equation (||) does not admit reg- 
ular self-similar solutions, so the numerically observed 
self-similarity can be only approximate. We identify the 
blowup proflle as the scale-evolving static solution (see 
Fig. 3) 
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We call this solution static because for any fixed A it is the 
time-independent solution of equation (|^) (this solution 
is perhaps better known as the YM instanton in four 
euclidean dimensions). Since the energy does not depend 
on A, these solutions are only marginally stable: when 
kicked they shrink or expand. In other words the blowup 
can be viewed as the adiabatic shrinking of the static 
solution. Numerical evidence suggests that the rate of 
blowup goes asymptotically to zero, that is ( ' = d/dt) 
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Although we are not able to explain this fact, we point 
out that it seems necessary for the consistency of the 
quasi-static character of blowup (||) . To see this, substi- 
tute (^ into (Q) to obtain 



(1 - X^ri^)W" + l + iXX- 2)?)if 

+ \w{i - W^) = 0. 
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It follows fro m (|lC| ) that the terms involving time deriva- 
tives of A in (|ll|) become asymptotically negligible and 
therefore in the leading order this equation has the same 
form as the right-hand side of (|^), which explains why 
the blowup profile has the shape of the static solution. 

In view of the limited resolution of our numerics and 
the lack of theory we are not in position to make con- 
jectures about the time dependence of A which would 
go beyond equation (|l^). Although a power law fit 
A ~ (T — t)^~^" with the anomalous exponent a « 0.1 
is quite accurate, we would not take this fact too seri- 
ously because we cannot rule out logarithmic corrections 
and, moreover, the exponent a exhibits weak dependence 
on initial data. 
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FIG. 3. Blowup in d = 4. The plot shows the late time 
evolution of large initial data (^4 = 0.5). The inner solution 
has the form of the scale-evolving static solution Ws{r/X{t)) 
with the scale X{t) going to zero slightly faster than linearly. 

It is instructive to compare the d = 4 and d — 5 
blowups from the standpoint of energy concentration. To 
this end, for solutions which blowup we define the energy 
at time t < T inside the past light cone of the singularity 
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where the coefficient c{d) = {d—l)vol{S'''^^) follows from 
integrating (j2|) over the angles and taking the trace. For 
d — 5, substituting (^ into ( p^ ) we obtain linij^T ^(i) = 
0, hence no energy gets concentrated into the singularity. 
In contrast, for c? = 4, substituting (||) into ( |l2|) and using 
( |l0| ) we obtain 
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thus the energy equal to the energy of the static solu- 
tion Ws concentrates at the singularity. Note that only 
the potential energy becomes concentrated; the kinetic 



energy tends asymptotically to zero. In fact, in our sim- 
ulations we can see the excess energy being slowly ra- 
diated away from the inner region as the blowup profile 
converges to the static solution. 

To summarize, our work provides numerical evidence 
that solutions of YM equations in four and five spatial di- 
mensions do form singularities from generic smooth large 
initial data. While the self-similar character of blowup 
in = 5 is well understood (at least from the numerical 
perspective) , the d = A case is more subtle and our anal- 
ysis leaves two important questions open, namely: what 
is the precise rate of blowup and what is the nature of 
the threshold for blowup? We plan to approach these 
issues by interpreting the expression (H) in terms of mo- 
tion along the one-dimensional moduli space of static so- 
lutions with the scale A playing the role of the collective 
coordinate. The most straightforward way of comput- 
ing the dynamics on the moduli space using the geodesic 
approximation is too naive in the present case: it pre- 
dicts that A changes linearly with time ||l^ which contra- 
dicts our numerics. We hope that more refined methods 
of dealing with collective coordinates, like the ones de- 
scribed in [p^ in the context of the nonlinear Schrodinger 
equation in two spatial dimensions, can be applied to our 
problem as well. In our opinion, the derivation of the cor- 
rect modulation equation for the scale A is the most im- 
portant next step towards understanding the dynamics 
of blowup for YM equations in four spatial dimensions; 
hopefully it would also shed light on the character of 
transition between blowup and dispersion. 
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